Redox Properties of Bacillus subtilis Ferredoxin:NADP+ Oxidoreductase: Potentiometric Characteristics and Reactions with Pro-Oxidant Xenobiotics

Bacillus subtilis ferredoxin:NADP+ oxidoreductase (BsFNR) is a thioredoxin reductase-type FNR whose redox properties and reactivity with nonphysiological electron acceptors have been scarcely characterized. On the basis of redox reactions with 3-acetylpyridine adenine dinucleotide phosphate, the two-electron reduction midpoint potential of the flavin adenine dinucleotide (FAD) cofactor was estimated to be −0.240 V. Photoreduction using 5-deazaflavin mononucleotide (5-deazaFMN) as a photosensitizer revealed that the difference in the redox potentials between the first and second single-electron transfer steps was 0.024 V. We examined the mechanisms of the reduction of several different groups of non-physiological electron acceptors catalyzed by BsFNR. The reactivity of quinones and aromatic N-oxides toward BsFNR increased when increasing their single-electron reduction midpoint redox potentials. The reactivity of nitroaromatic compounds was lower due to their lower electron self-exchange rate, but it exhibited the same trend. A mixed single- and two-electron reduction reaction was characteristic of quinones, whereas reactions involving nitroaromatics proceeded exclusively via the one-electron reduction reaction. The oxidation of FADH• to FAD is the rate-limiting step during the oxidation of fully reduced FAD. The calculated electron transfer distances in the reaction with nitroaromatics were close to those of other FNRs including the plant-type enzymes, thus demonstrating their similar active site accessibility to low-molecular-weight oxidants despite the fundamental differences in their structures.


Introduction
Ferredoxin:NADP + oxidoreductases (FNRs, EC 1.18.1.2) are ubiquitous flavoenzymes found in all domains of life and belong to the class of dehydrogenase-electrontransferases.In oxygenic photosynthesis, monomeric plant-type FNR with a molecular mass of 35-36 kDa transfers a redox equivalent from reduced iron-sulfur protein ferredoxin (Fd) to NADP + , providing NADPH for CO 2 assimilation [1][2][3][4][5].In non-photosynthetic processes, FNR catalyzes the reduction of Fd at the expense of NADPH supplying reduced Fd for nitrate assimilation (roots), biosynthesis of isoprenoids (the malaria parasite Plasmodium falciparum), and so on [5][6][7].Under iron-limiting conditions, the FMN-containing low-potential electron carriers flavodoxins (Flds) function as the physiological redox partners of FNR in place of Fd [8,9].However, stronger and more specific interactions occur between FNR and Fd in comparison to those occurring with Fld [10].
FNR from the Gram-positive bacterium Bacillus subtilis (BsFNR) belongs to a novel group of FNRs classified as NADPH-thioredoxin reductase (TrxR)-type FNR based on their high sequence similarity with bacterial TrxRs [3,11].A feature that distinguishes B. subtilis, as well as other representatives of the TrxR-type FNR, FNRs from Rhodopseudomonas palustris, and Chlorobaculum tepidum, from plant-type FNR is their homodimeric structures (Figure 1) [12,13].Another distinguishing feature is that the FAD-binding domain comprises two discontinuous segments separated by the segment of the NADPH-binding domain (Figure 1), while in plant-type FNRs, each domain is composed of a single peptide segment [3].In BsFNR, the FAD-binding domain is composed of residues 1-124 and 250-332, and the NADP(H)-binding domain is composed of residues 125-249.The FADbinding and NADP(H)-binding domains are connected by two flexible hinge regions that enable the rotation of the domains relative to one another.This is plausible since, in the case of BsFNR, the distance between bound NADP + and the isoalloxazine ring of FAD is ca.15 Å in the crystal structure, which is too distant for an efficient hydride transfer.The crystal structure of BsFNR shows two aromatic amino acids, Tyr50 and His324, located over the si and re faces of the isoalloxazine ring of FAD, respectively.These amino acids are almost parallel to the isoalloxazine ring, with the corresponding distance being 3.3-3.6Å, and form a π-π stacked structure [13].
group of FNRs classified as NADPH-thioredoxin reductase (TrxR)-type FNR based on their high sequence similarity with bacterial TrxRs [3,11].A feature that distinguishes B. subtilis, as well as other representatives of the TrxR-type FNR, FNRs from Rhodopseudomonas palustris, and Chlorobaculum tepidum, from plant-type FNR is their homodimeric structures (Figure 1) [12,13].Another distinguishing feature is that the FADbinding domain comprises two discontinuous segments separated by the segment of the NADPH-binding domain (Figure 1), while in plant-type FNRs, each domain is composed of a single peptide segment [3].In BsFNR, the FAD-binding domain is composed of residues 1-124 and 250-332, and the NADP(H)-binding domain is composed of residues 125-249.The FAD-binding and NADP(H)-binding domains are connected by two flexible hinge regions that enable the rotation of the domains relative to one another.This is plausible since, in the case of BsFNR, the distance between bound NADP + and the isoalloxazine ring of FAD is ca.15 Å in the crystal structure, which is too distant for an efficient hydride transfer.The crystal structure of BsFNR shows two aromatic amino acids, Tyr50 and His324, located over the si and re faces of the isoalloxazine ring of FAD, respectively.These amino acids are almost parallel to the isoalloxazine ring, with the corresponding distance being 3.3-3.6Å, and form a π-π stacked structure [13].
Figure 1.The 3D structure of a BsFNR dimer (PDB id: 3LZX).FAD-binding domain is colored yellow, while the NADP(H)-binding domain is colored red.The box denotes amino acids comprising each domain in a monomer.The bound FAD is shown in purple.The molecular graphics were made using UCSF ChimeraX [14].
The mechanism of the catalysis of BsFNR with its physiological substrates has been partially characterized.At pH 7.0 and 10 °C, NADPH reduces the enzyme, with k ≥ 500 s −1 , while its reoxidation with NADP + occurs much more slowly and partially [15].The reactions are accompanied by the transient formation of charge transfer complexes of NADP + /H with reduced/oxidized FAD, respectively, which absorb at longer wavelengths.The maximal steady-state turnover rates of BsFNR using the non-physiological electron acceptor ferricyanide or Fe4S4-containing B. subtilis ferredoxin (BsFd) are close to 1000 s −1 and 50-100 s −1 based on single-electron transfer, respectively [16,17].In the reduction of BsFNR with reduced BsFd, the rate-limiting step is the reduction of the FAD semiquinone (FADH • ) intermediate, with k > 50 s −1 [15].A possible physiological function of BsFNR is the reduction of BsFd and flavodoxins, which further act as electron donors for cytochrome P450 BioI [18,19], lipid desaturase [20], and NO synthase, which, in this case, consists only of the oxygenase domain [21].The mechanism of the catalysis of BsFNR with its physiological substrates has been partially characterized.At pH 7.0 and 10 • C, NADPH reduces the enzyme, with k ≥ 500 s −1 , while its reoxidation with NADP + occurs much more slowly and partially [15].The reactions are accompanied by the transient formation of charge transfer complexes of NADP + /H with reduced/oxidized FAD, respectively, which absorb at longer wavelengths.The maximal steady-state turnover rates of BsFNR using the non-physiological electron acceptor ferricyanide or Fe 4 S 4 -containing B. subtilis ferredoxin (BsFd) are close to 1000 s −1 and 50-100 s −1 based on single-electron transfer, respectively [16,17].In the reduction of BsFNR with reduced BsFd, the rate-limiting step is the reduction of the FAD semiquinone (FADH • ) intermediate, with k > 50 s −1 [15].A possible physiological function of BsFNR is the reduction of BsFd and flavodoxins, which further act as electron donors for cytochrome P450 BioI [18,19], lipid desaturase [20], and NO synthase, which, in this case, consists only of the oxygenase domain [21].
There are other, virtually untouched aspects of BsFNR redox reactions.In particular, B. subtilis, as a soil bacterium, is naturally exposed to redox-active humic compounds, quinones generated by various fungi and bacteria, and nitroaromatic industrial com-pounds [22][23][24][25].Like other FNRs, BsFNR is a potential target of these compounds, and it can stimulate their toxic effects by reducing them to free radicals that cause oxidative stress, or, alternatively, it can participate in their reductive detoxification [25,26].Another aspect is the possible action of BsFNR as a target of redox-cycling drugs, analogous to what was suggested for the malaria parasite Plasmodium falciparum FNR [27].Although B. subtilis is generally harmless to humans, separate cases of it being pathogenic have been reported [28,29].On the other hand, TrxR-type FNRs have been found in most firmicutes, including pathogens, and, for example, BsFNR has 70% amino acid sequence homology with TrxR-like FNR from pathogenic Listeria monocytogenes [30].For these reasons, we employed steady-and pre-steady-state kinetics methods to determine the structure-activity relationships of BsFNR in reactions with a series of redox-cycling quinones, nitroaromatic compounds, and aromatic N-oxides.For the quantitative analysis of the obtained data and comparison with other FNRs, the redox potentials of BsFNR were also determined.

The Determination of Redox Potentials of BsFNR
The potentiometric properties of TrxR-type FNRs are largely unexplored, with the exception of the previously determined standard redox potential (E 0 7 , the potential of the E-FAD/E-FADH − redox couple) and corresponding single-electron transfer potentials (E-FAD/E-FADH • and E-FADH • /E-FADH redox couples) of R. palustris FNR [31].The E 0 7 of flavoenzymes can be determined using the Haldane relationship, according to which the ratio of the bimolecular rate constants of forward and reverse reactions yields the equilibrium constant of the reaction (K) [31,32].This, in turn, is related to the difference in the standard redox potential of the reactants (∆E 0 (V) = 0.0295 × log K for a two-electron transfer).Because the estimation of the kinetic parameters of the reverse reaction, the reduction of NADP + by BsFNR, is complicated by BsFd substrate inhibition [15], we examined the enzyme reactions using the NADP(H) analogue 3-acetylpyridine adenine dinucleotide phosphate (APADP(H), E 0 7 = −0.258V [33]).In this case, K is expressed as the ratio of bimolecular rate constants (k cat /K m ) of APADPH oxidation and APADP + reduction.In the forward reaction, APADPH (50-500 µM) was generated in situ using the glucose-6phosphate/glucose-6-phosphate dehydrogenase reduction system, and the reaction rate was monitored following the reduction of 1.0 mM of ferricyanide.The reaction rates did not depend on ferricyanide concentration in the range of 0.25-1.0mM.The reaction was characterized by k app cat = 105.0± 3.2 s −1 and k cat /K m = 1.5 ± 0.2 × 10 6 M −1 s −1 (on a two-electron basis).The enzymatic reduction of APADP + (50-500 µM) by 200 µM of NADPH proceeded with k app cat = 35.9± 1.0 s −1 and k cat /K m = 3.6 ± 0.4 × 10 5 M −1 s −1 .Thus, according to calculations K = 4.1 ± 0.7, and E 0 7 for the enzyme was −0.240 ± 0.002 V.During the photoreduction of BsFNR with 5-deazaFMN and EDTA as a photosensitizer, a neutral (blue) FAD semiquinone (FADH • ) transiently formed, as evidenced by the characteristic broad absorption band in the 550-650 nm range (Figure 2).The extinction coefficient ε 600 for BsFNR FAD is not known; however, the amount of E-FADH • can be assessed from the A 457 vs.A 600 plot (Figure 2 inset) using the value ε 600 = 3.6 mM −1 cm −1 for the Y50G BsFNR mutant [34].In this case, the maximal amount of semiquinone stabilized was 44%.The separation between the two single-electron-transfer potentials (∆E 1  7 ) can be calculated based on the semiquinone formation constant K s (Equations ( 1) and (2)): where [E − FADH • ] max is the maximal concentration of semiquinone and [E − FAD] tot is the total concentration of enzyme.The calculations yield K s = 2.572 and ∆E 1 7 = 0.024 V, and, subsequently, E E−FAD/E−FADH •

Steady-State Kinetics and Oxidant Substrate Specificity of BsFNR
Earlier studies performed on FNRs from other sources indicate that juglone (5hydroxy-1,4-naphthoquinone) and its derivatives are potent nonphysiological oxidants of these enzymes [31,35].A series of parallel lines were obtained in double reciprocal plots upon varying the concentration of NADPH with a constant concentration of juglone (Figure 3) and vice versa.This shows that the BsFNR-catalyzed quinone reductase reaction proceeds via a "ping-pong" mechanism.

Steady-State Kinetics and Oxidant Substrate Specificity of BsFNR
Earlier studies performed on FNRs from other sources indicate that juglone (5-hydroxy-1,4-naphthoquinone) and its derivatives are potent nonphysiological oxidants of these enzymes [31,35].A series of parallel lines were obtained in double reciprocal plots upon varying the concentration of NADPH with a constant concentration of juglone (Figure 3) and vice versa.This shows that the BsFNR-catalyzed quinone reductase reaction proceeds via a "ping-pong" mechanism.

Steady-State Kinetics and Oxidant Substrate Specificity of BsFNR
Earlier studies performed on FNRs from other sources indicate tha hydroxy-1,4-naphthoquinone) and its derivatives are potent nonphysiologic these enzymes [31,35].A series of parallel lines were obtained in double re upon varying the concentration of NADPH with a constant concentratio (Figure 3) and vice versa.This shows that the BsFNR-catalyzed quinone redu proceeds via a "ping-pong" mechanism.As calculated according to Equation (9) (presented in the Materials and Methods Section), the k cat value for juglone reduction at an infinite NADPH concentration is equal to 360.2 ± 9.8 s −1 , and the values of the bimolecular rate constants (k cat /K m ) for juglone and NADPH are equal to 4.35 ± 0.23 × 10 6 M −1 s −1 and 1.00 ± 0.07 × 10 7 M −1 s −1 , respectively.
The quinone reductase reaction of BsFNR is inhibited by the reaction product NADP + .At a fixed oxidant concentration, NADP + acts as a competitive inhibitor towards NADPH, increasing the slopes of the Lineweaver-Burk plots but not affecting the maximal rate of reaction (Figure 4A).Its inhibition constant (K i ), calculated according to Equation (10) (Materials and Methods), is equal to 62.6 ± 10.5 µM.In turn, at a fixed NADPH concentration, NADP + acts as an uncompetitive inhibitor towards the oxidant, decreasing the maximal reaction rate but not increasing the slopes of the Lineweaver-Burk plots (Figure 4B).Its uncompetitive K i , calculated according to Equation ( 11) (shown in the Materials and Methods Section), is equal to 261 ± 27 µM.
Int. J. Mol.Sci.2024, 25, x 5 of 17 As calculated according to Equation ( 9) (presented in the Materials and Methods Section), the kcat value for juglone reduction at an infinite NADPH concentration is equal to 360.2 ± 9.8 s −1 , and the values of the bimolecular rate constants (k cat K m ⁄ ) for juglone and NADPH are equal to 4.35 ± 0.23 × 10 6 M −1 s −1 and 1.00 ± 0.07 × 10 7 M −1 s −1 , respectively.The quinone reductase reaction of BsFNR is inhibited by the reaction product NADP + .At a fixed oxidant concentration, NADP + acts as a competitive inhibitor towards NADPH, increasing the slopes of the Lineweaver-Burk plots but not affecting the maximal rate of reaction (Figure 4A).Its inhibition constant (Ki), calculated according to Equation (10) (Materials and Methods), is equal to 62.6 ± 10.5 µM.In turn, at a fixed NADPH concentration, NADP + acts as an uncompetitive inhibitor towards the oxidant, decreasing the maximal reaction rate but not increasing the slopes of the Lineweaver-Burk plots (Figure 4B).Its uncompetitive Ki, calculated according to Equation (11) (shown in the Materials and Methods Section), is equal to 261 ± 27 µM.A series of quinones (Q), aromatic nitrocompounds (ArNO2), and aromatic N-oxides (ArN→O) with single-electron reduction potentials ( E ) from 0.09 to −0.494 V were studied to determine the specificity of BsFNR towards oxidizing substrates.The latter group of compounds, derivatives of 3-amino-1,2,4-benzotriazine-1,4-dioxide (tirapazamine), were studied in connection with their antimicrobial activity [36].The single-electron acceptors ferricyanide, benzyl viologen, and FeEDTA − were also studied.Among the studied forms of ArNO2, there are explosives, such as tetryl and 2,4,6trinitrotoluene; antibacterial agents, like nitrofurantoin and nifuroxime; and the anticancer agent CB-1954.The apparent maximal reduction rate constants k cat app of the electron acceptors at 100 µM of NADPH and their respective k cat K m ⁄ values are given in Table 1.A series of quinones (Q), aromatic nitrocompounds (ArNO 2 ), and aromatic N-oxides (ArN→O) with single-electron reduction potentials (E 1  7 ) from 0.09 to −0.494 V were studied to determine the specificity of BsFNR towards oxidizing substrates.The latter group of compounds, derivatives of 3-amino-1,2,4-benzotriazine-1,4-dioxide (tirapazamine), were studied in connection with their antimicrobial activity [36].The single-electron acceptors ferricyanide, benzyl viologen, and FeEDTA − were also studied.Among the studied forms of ArNO 2 , there are explosives, such as tetryl and 2,4,6-trinitrotoluene; antibacterial agents, like nitrofurantoin and nifuroxime; and the anticancer agent CB-1954.The apparent maximal reduction rate constants k app cat of the electron acceptors at 100 µM of NADPH and their respective k cat /K m values are given in Table 1.
The log k cat /K m values for ArNO 2 exhibit a linear dependence on their E 1 7 values (Table 1, Figure 5), with tetryl (16) being the most potent electron acceptor in that group.On the other hand, the log k cat /K m values for quinones with a similar E 1 7 are about an order of magnitude higher.Moreover, quinones and aromatic N-oxides together exhibit a distinct parabolic dependence on their E 1 7 values (Figure 5).One should also note the single-electron acceptor benzyl viologen (40), with its steady-state kinetic constants being close to those of quinones with similar E 1 7 values.We found that BsFNR reduces Q and ArNO 2 to their radicals.Typically, the singleelectron flux for quinone reduction by NAD(P)H-oxidizing enzymes is defined as the ratio of the rate of 1,4-benzoquinone-mediated reduction of the added cytochrome c to the double rate of 1,4-benzoquinone-mediated NAD(P)H enzymatic oxidation at pH < 7.2 [41].Such an approach is based on the fast reduction of cytochrome c by the 1,4-benzosemiquinone (k ≈ 10 6 M −1 s −1 ) and its slow reduction by the hydroquinone form.In the case of the BsFNR-catalyzed reduction of 20-100 µM of 1,4-benzoquinone by 100 µM of NADPH, the single-electron flux was equal to 17%.On the other hand, during the reduction of 9,10- The log k cat K m ⁄ values for ArNO2 exhibit a linear dependence on their E values (Table 1, Figure 5), with tetryl (16) being the most potent electron acceptor in that group.On the other hand, the log k cat K m ⁄ values for quinones with a similar E are about an order of magnitude higher.Moreover, quinones and aromatic N-oxides together exhibit a distinct parabolic dependence on their E values (Figure 5).One should also note the single-electron acceptor benzyl viologen (40), with its steady-state kinetic constants being close to those of quinones with similar E values.1.
We found that BsFNR reduces Q and ArNO2 to their radicals.Typically, the singleelectron flux for quinone reduction by NAD(P)H-oxidizing enzymes is defined as the ratio of the rate of 1,4-benzoquinone-mediated reduction of the added cytochrome c to the double rate of 1,4-benzoquinone-mediated NAD(P)H enzymatic oxidation at pH < 7.2 [41].Such an approach is based on the fast reduction of cytochrome c by the 1,4benzosemiquinone (k ≈ 10 6 M −1 s −1 ) and its slow reduction by the hydroquinone form.In the case of the BsFNR-catalyzed reduction of 20-100 µM of 1,4-benzoquinone by 100 µM of NADPH, the single-electron flux was equal to 17%.On the other hand, during the reduction of 9,10-phenanthrenequinone, the calculated single-electron flux was equal to 70%, with superoxide dismutase (SOD) (100 U/mL) inhibiting cytochrome c reduction by 43%.The assessment of the single-electron flux for the reduction of ArNO2 is also based on the ArNO2 •− -mediated reduction of added cytochrome c.In the presence of 50 µM of NADPH and 100 µM of TNT or p-nitrobenzaldehyde, the rate of BsFNR-catalyzed reduction of 50 µM of cytochrome c was equal to 180-200% of the NADPH oxidation rate, corresponding to 90-100% single-electron flux.SOD inhibited the reactions by 10-30%, reflecting the rapid reoxidation of the radicals with O2 and the participation of the superoxide in the reduction of cytochrome c.
In order to characterize the role of electrostatic interactions in the reduction of oxidants, experiments were performed, with varying ionic strengths, using negatively charged oxidant ferricyanide, positively charged benzyl viologen, and neutral 1,4naphthoquinone (Figure 6).The log k cat K m ⁄ for ferricyanide exhibits a parabolic character with a peak value at I ≈ 0.25 M. The reactivity of uncharged 1,4-napthoquinone decreases Figure 5. Dependence of the reactivity of quinones, nitroaromatic compounds, aromatic N-oxides, and benzyl viologen on their single-electron reduction midpoint potentials (log 10 scale).Numbers and reduction potentials of the compounds are given in Table 1.
In order to characterize the role of electrostatic interactions in the reduction of oxidants, experiments were performed, with varying ionic strengths, using negatively charged oxidant ferricyanide, positively charged benzyl viologen, and neutral 1,4-naphthoquinone (Figure 6).The log k cat /K m for ferricyanide exhibits a parabolic character with a peak value at I ≈ 0.25 M. The reactivity of uncharged 1,4-napthoquinone decreases slightly upon increasing ionic strength.In general, the log k cat /K m of benzyl viologen decreased upon increasing the ionic strength of the buffer solution (Figure 6).slightly upon increasing ionic strength.In general, the log k cat K m ⁄ of benzyl viologen decreased upon increasing the ionic strength of the buffer solution (Figure 6).

Pre-steady-State Kinetics of BsFNR Oxidation under Multiple Turnover Conditions
Some insight into the reoxidation mechanism of BsFNR can be gleaned from the spectral changes of the enzyme-bound FAD during the multiple turnovers under aerobic conditions in the presence of NADPH and tetramethyl-1,4-benzoquinone (duroquinone,

Pre-Steady-State Kinetics of BsFNR Oxidation under Multiple Turnover Conditions
Some insight into the reoxidation mechanism of BsFNR can be gleaned from the spectral changes of the enzyme-bound FAD during the multiple turnovers under aerobic conditions in the presence of NADPH and tetramethyl-1,4-benzoquinone (duroquinone, DQ).DQ does not absorb above 460 nm, and its semiquinone form is rapidly reoxidized by oxygen [37].The control experiment performed in the absence of DQ showed an initial fast phase of FAD reduction by NADPH, which was then followed by a much slower reoxidation by O 2 seen at 460 nm and a transient increase in absorbance at 600 nm (Figure 7A), at the same time scale.The addition of DQ enhanced the reoxidation rate by up to two orders of magnitude (Figure 7B).

Pre-steady-State Kinetics of BsFNR Oxidation under Multiple Turnover Conditions
Some insight into the reoxidation mechanism of BsFNR can be gleaned from the spectral changes of the enzyme-bound FAD during the multiple turnovers under aerobic conditions in the presence of NADPH and tetramethyl-1,4-benzoquinone (duroquinone, DQ).DQ does not absorb above 460 nm, and its semiquinone form is rapidly reoxidized by oxygen [37].The control experiment performed in the absence of DQ showed an initial fast phase of FAD reduction by NADPH, which was then followed by a much slower reoxidation by O2 seen at 460 nm and a transient increase in absorbance at 600 nm (Figure 7A), at the same time scale.The addition of DQ enhanced the reoxidation rate by up to two orders of magnitude (Figure 7B).The reaction rate constants were calculated according to Equation (3) [42]: where [NADPH]0 is the initial concentration of NADPH, [Ered]max is the maximal concentration of the reduced enzyme formed during the turnover, and t1/2(off) is the time interval between the formation of the half-maximum amount of Ered and its decay to the half-maximum value.According to the absorbance difference between the oxidized and The reaction rate constants were calculated according to Equation (3) [42]: where [NADPH] 0 is the initial concentration of NADPH, [E red ] max is the maximal concentration of the reduced enzyme formed during the turnover, and t 1/2(off) is the time interval between the formation of the half-maximum amount of E red and its decay to the half-maximum value.According to the absorbance difference between the oxidized and reduced enzyme forms, that is, ∆ε 460 = 9.02 mM −1 cm −1 [15], the maximal extent of enzyme reduction is 91% (Figure 7A).This yields k ox = 0.22 s −1 for the reoxidation of BsFNR by O 2 , a value that is similar to the enzyme NADPH oxidase activity, 0.17 s −1 .Analogous calculations at various DQ concentrations (Figure 8) yield the value of k ox(max .) = 48.3± 6.7 s −1 (Figure 8 inset), which may reflect the maximum rate of the oxidative half-reaction.
Int. J. Mol.Sci.2024, 25, x 9 of 17 reduced enzyme forms, that is, Δε460 = 9.02 mM −1 cm −1 [15], the maximal extent of enzyme reduction is 91% (Figure 7A).This yields kox = 0.22 s −1 for the reoxidation of BsFNR by O2, a value that is similar to the enzyme NADPH oxidase activity, 0.17 s −1 .Analogous calculations at various DQ concentrations (Figure 8) yield the value of kox(max.)= 48.3± 6.7 s −1 (Figure 8 inset), which may reflect the maximum rate of the oxidative half-reaction.Measurements taken at different wavelengths showed that the intermediate product Measurements taken at different wavelengths showed that the intermediate product of reoxidation by DQ has a flat absorption maximum between 600 and 700 nm (Figure 9).The decay of the reaction intermediate begins at about 500 ms after the start of the reaction.

Figure 8.
The kinetics of BsFNR (5 µM) reduction and reoxidation under multiple turnover conditions with a varying concentration of duroquinone in the presence of 50 µM of NADPH followed at 460 nm (solid curves) and 600 nm (dashed curves): 1-250 µM of duroquinone, 2-166 µM of duroquinone, 3-111 µM of duroquinone, 4-74 µM of duroquinone, 5-49 µM of duroquinone, and 6-33 µM of duroquinone (all concentrations were reported after mixing).The inset shows the Lineweaver-Burk plot for the dependence of the apparent first-order rate constant of the reoxidation of reduced FAD on the concentration of duroquinone.
Measurements taken at different wavelengths showed that the intermediate product of reoxidation by DQ has a flat absorption maximum between 600 and 700 nm (Figure 9).The decay of the reaction intermediate begins at about 500 ms after the start of the reaction.

Discussion
Since TrxR-type FNRs are relatively poorly studied enzymes, our studies focused on revealing similarities and differences between the redox properties of BsFNR and other groups of FNRs.For the first time, we determined the standard redox potential of BsFNR, i.e., −0.240 V, which is significantly more positive than the E 0 7 of FNRs from other sources, e.g., −0.280 V (P.falciparum [43]) or −0.342 V (spinach [44]).Most importantly, however, it is sufficiently more positive than the potential of its TrxR-type homologue, FNR from R. palustris, with a value of −0.276 V [31].The amino acids surrounding the isoalloxazine ring of FAD in both enzymes are largely conserved, except for the change of Tyr328 to His324 and Thr329 to Ser325 in BsFNR (Figure 10) [12,13].The possible effect of the latter substitution on BsFNR E 0 7 is not clear, but the effect of the presence of His324 is similar to the case for flavodoxin from Desulfovibrio vulgaris , where Tyr98His substitution increased the E 0 7 of FMN by 0.07 V [45].It is believed that the presence of the imidazole group of histidine stabilizes the anionic form of the reduced flavin, i.e., makes its oxidation more difficult.This elevated E 0 7 value agrees with previous data on the inefficient reoxidation of BsFNR by NADP + [15] and the possible role of this enzyme in Fd/Fld reduction at the expense of NADPH.On the other hand, the FAD semiquinone stability of BsFNR, 44% at equilibrium (Figure 2), is slightly but not drastically higher than that of R. palustris FNR, i.e., 26.5% [31].The stability of FADH • may be due to the presence of Tyr50, Asp57, Ile296, and Thr326, the last of which forms a H-bond with N5 of isoalloxazine, in the active site of BsFNR, corresponding to the Tyr49, Asp56, Ile 300, and Thr330 of R. palustris FNR [46].The semiquinone-stabilizing effect of the aforementioned amino acid residues has been discussed in our previous work [31].
of BsFNR, 44% at equilibrium (Figure 2), is slightly but not drastically higher than that of R. palustris FNR, i.e., 26.5% [31].The stability of FADH • may be due to the presence of Tyr50, Asp57, Ile296, and Thr326, the last of which forms a H-bond with N5 of isoalloxazine, in the active site of BsFNR, corresponding to the Tyr49, Asp56, Ile 300, and Thr330 of R. palustris FNR [46].The semiquinone-stabilizing effect of the aforementioned amino acid residues has been discussed in our previous work [31].Due to the possible involvement of BsFNR in reactions with redox-active agents in the soil, further attention was paid to the analysis of reactions of model redox-active compounds of different groups.Here, we found a series of features common to a wide range of FNRs: (i) the "ping-pong" mechanism of reaction and the way of NADP + inhibition (Figures 3 and 4) are similar to those of other FNRs [31,35].This result shows that the reductive and oxidative half-reactions occur separately but at the same binding site.The considerable variation in the k cat app values of the reactions with different oxidants (Table 1) suggest that oxidative half-reaction is the rate-limiting step.The competitive inhibition of NADP + with respect to NADPH and the uncompetitive inhibition with respect to the oxidant is explained by the specific case of the "ping-pong" mechanism, wherein NADP + strongly competes with NADPH for binding to the oxidized enzyme but binds poorly to the reduced form of the enzyme and/or does not reoxidize it [48].(ii) Another feature common to various classes of FNRs is the dependence of log k cat K m ⁄ of various classes of oxidants on their E , i.e., the absence of a pronounced impact of their structural features on reactivity (Figure 5) [31,35].This may correspond to the "outer sphere" single-electron transfer model [49].In this case, the rate constant for the electron Due to the possible involvement of BsFNR in reactions with redox-active agents in the soil, further attention was paid to the analysis of reactions of model redox-active compounds of different groups.Here, we found a series of features common to a wide range of FNRs: (i) the "ping-pong" mechanism of reaction and the way of NADP + inhibition (Figures 3 and 4) are similar to those of other FNRs [31,35].This result shows that the reductive and oxidative half-reactions occur separately but at the same binding site.The considerable variation in the k app cat values of the reactions with different oxidants (Table 1) suggest that oxidative half-reaction is the rate-limiting step.The competitive inhibition of NADP + with respect to NADPH and the uncompetitive inhibition with respect to the oxidant is explained by the specific case of the "ping-pong" mechanism, wherein NADP + strongly competes with NADPH for binding to the oxidized enzyme but binds poorly to the reduced form of the enzyme and/or does not reoxidize it [48].(ii) Another feature common to various classes of FNRs is the dependence of log k cat /K m of various classes of oxidants on their E 1 7 , i.e., the absence of a pronounced impact of their structural features on reactivity (Figure 5) [31,35].This may correspond to the "outer sphere" single-electron transfer model [49].In this case, the rate constant for the electron transfer between reactants (k 12 ) depends on their electron self-exchange rate constants (k 11 and k 22 ) and the equilibrium constant of reaction (K): and f is expressed as where Z is the frequency factor, 10 11 M −1 s −1 .According to Equations ( 5) and (6) for the reaction between an electron donor and a series of homologous electron acceptors with similar k 22 values, log k 12 will exhibit a parabolic (or linear, in the case of low exothermicity) dependence on ∆E 1 .We observed this experimentally: the k 22 of ArNO 2 , ≈ 10 6 M −1 s −1 , is 100 times lower than that of Q and ArN→O, ≈ 10 8 M −1 s −1 [40,50], leading to a 10-fold lower reactivity of ArNO 2 when compared to quinones of similar E 1 7 values (Figure 5).(iii) The transient formation of 500-700 nm absorbing species during enzyme reoxidation is typical for other FNRs [31,35] and consistent with a two-step reoxidation scheme, namely, FADH − →FADH • →FAD, with FADH • oxidation being the rate-limiting step.This means that the obtained k cat /K m of the compounds (Table 1) reflect FADH • oxidation.The small absorption increase above 700 nm is uncharacteristic of FADH • and indicative of a parallel formation of other reaction intermediates, e.g., the complexes of FADH • with NADP(H) observed in adrenodoxin reductase [51][52][53][54].
In parallel, some differences from previously studied FNRs were also observed.The experiments conducted at varied ionic strengths may characterize the region of BsFNR that interacts with charged oxidants.The data on Anabaena PCC7119 and P. falciparum FNRs demonstrate the bell-shaped dependences of log k cat /K m on ionic strength irrespective of the charge of the oxidants [35].This can be attributed to an interaction with both the negatively charged glutamate and positively charged lysine residues, which are located close to the dimethylbenzene part of the isoalloxazine ring [55][56][57].However, the data in Figure 6 show that the oxidants interact with the negatively charged amino acid residue(s) of BsFNR.In this case, the potential candidates are Asp285 and, possibly, Asp57 near the isoalloxazine ring of FAD [13].Most likely, the parallel conformational changes in the protein occur at high ionic strength, decreasing the reactivity of the uncharged 1,4-naphthoquinone (Figure 6).The positively charged fragment of the BsFNR FAD domain, Lys264, Lys317, andArg319 [13], which possibly electrostatically interacts with the negatively charged BsFd [18], may not interact with non-physiological electron acceptors because it is too far from the isoalloxazine ring [15].However, the most unexpected difference from other FNRs was the mixed one-and two-electron reduction of quinones by BsFNR.Both Anabaena PCC7119 and P. falciparum FNRs reduce both quinones and ArNO 2 only in a single-electron fashion [35,58].The fact that the Glu301Ala mutant of Anabaena FNR reduced 1,4-benzoquinone in a 50% single-electron manner can be explained by its significant destabilization of FADH • , 8% at equilibrium [58].However, BsFNR reduces 1,4-benzoquinone with only 17% of single-electron flux, although the FADH • stabilization is 44% (Figure 2).Similarly, another TrxR-type R. palustris FNR with 26.5% FADH • stabilization reduces 1,4-benzoquinone with 54% single-electron flux [31].Since this regularity was observed in the case of two TrxR-type FNRs, it is possible that single-and two-electron reduction flux are not determined by FADH • stability alone but also by the specificity of the isoalloxazine surroundings.One of the operative factors could be the unusual π-π stacking between the isoalloxazine ring and the aromatic Tyr50 and His324 groups [16,34].This feature, including its possible influence on the degradation of xenobiotics, is the subject of our further research.
Finally, the influence of BsFNR structural peculiarities in reactions with redox-active xenobiotics can be quantitatively described.According to Mauk et al., the electron transfer distance in metalloproteins reacting with inorganic complexes under infinite ionic strength (i.e., no electrostatic effects) can be related to the metalloproteins' electron self-exchange rate constant (k 11 ) (Equation ( 7)) [59]: R p () = 6.3 − 0.35ln k 11 (7) This approach has been used to estimate the R p values in the single-electron oxidation of various flavoenzyme dehydrogenase-electrontransferases [31,35].However, the values thus obtained should be regarded cautiously due to the possibility of electron transfer distances being overestimated since, in the case of dehydrogenase-electrontransferases, the dimethylbenzene part of the isoalloxazine ring is partly exposed to the solvent.Thus, these data can only be used to compare the "intrinsic" reactivity of different FNRs.
Since there is some uncertainty regarding the mechanism of quinone reduction by BsFNR, we will only analyze the case of ArNO 2 reduction.At ∆E , the potential of E-FADH • /E-FAD redox couple, and considering that the k 22 of ArNO 2 is equal to 10 6 M −1 s −1 [50], the value of log k 11 calculated according to the data in Figure 5 is equal to 3.11 ± 0.15.According to Equation (7), this yields R p = 3.8 ± 0.1 Å.For comparison, R p values calculated in an analogous way for nitroreductase reactions are equal to 4.4 Å (Anabaena FNR), 4.9-5.6Å (P.falciparum FNR), and 5.4 ± 0.2 Å (R. palustris FNR) [31,35].This shows that these R p s are close to each other, including for the plant-type FNRs.Complementing the R. palustris FNR data [31], this shows that the active sites of TrxR-type FNRs are not characterized by decreased accessibility to low-molecular-weight oxidants, which may be caused by the shielding of the isoalloxazine ring by the flexible C-terminal region, or by the domain movement of the protein.
In addition to providing fundamental information on the structure and catalysis of this class of FNRs, these data may also have some practical implications, e.g., with respect to using the FNRs of this group in bioelectrocatalysis [60].In this case, the environment of the redox cofactor and the distance from the surface of the protein usually determine the rate of the corresponding electrochemical reactions.

Steady-State Kinetics
The steady-state kinetics experiments were performed using a Cary60 UV/Vis spectrophotometer (Agilent Technologies, Santa Clara, CA, USA).All experiments were performed in 0.02 M Hepes/NaOH + 1 mM EDTA buffer solution, with a pH of 7.0, at 25 • C. The kinetic data were fitted to the Michaelis-Menten equation in Mathematica (Wolfram Research, Inc., Mathematica, Version 14.0, Champaign, IL, USA (2024)) (Equation ( 8)) to yield the steady-state parameters of the reactions, namely, the catalytic constants k app cat and bimolecular reaction rate constants (or catalytic efficiency constants) k cat /K m of the oxidants under a fixed concentration of NADPH: where v is the reaction rate, [E] is the concentration of BsFNR, [S] is the concentration of the oxidant, and k app cat represents the number of molecules of NADPH oxidized by a single native molecule of the enzyme per second at saturated concentrations of both substrates.The fitted parameters are equal to the reciprocal intercepts and slopes of Lineweaver-Burk plots, [E]/v vs. 1/[S], respectively.The concentrations of the enzyme used in these experiments were 5-50 nM.The kinetic parameters of the steady-state reactions according to the "ping-pong" mechanism were calculated according to Equation (9): The competitive inhibition constant K i of NADP + vs. NADPH was calculated according to Equation (10): The uncompetitive inhibition constant K i of NADP + vs. electron acceptor (Q) was calculated according to Equation (11): The rates of enzymatic NADPH oxidation in the presence of quinones, nitroaromatics, aromatic N-oxides, or single-electron acceptors were determined using the value ∆ε 340 = 6.2 mM −1 cm −1 , and they were corrected for the intrinsic NADPH-oxidase activity of BsFNR (0.17 s −1 ) and/or nonenzymatic NADPH oxidation by high-potential quinones.When 50 µM of cytochrome c was added to the reaction mixture, its quinone-and nitroaromatic-mediated reduction was assessed according to ∆ε 550 = 20 mM −1 cm −1 .The ferricyanide reduction rate was measured according to ∆ε 420 = 1.03 mM −1 cm −1 .The rate of BsFNR-catalyzed reduction of APADP + by NADPH was determined according to ∆ε 363 = 5.6 mM −1 cm −1 [33].APADPH was prepared in situ by reducing APADP + with 10 mM of glucose 6-phosphate and 0.01 mg/mL of glucose 6-phosphate dehydrogenase, and its concentration was determined according to ε 365 = 7.8 mM −1 cm −1 [33].NaCl was used to vary the ionic strength of the buffer solution.The stock solutions of organic compounds were prepared in DMSO; the final concentration of DMSO in reaction mixtures was 1% (v/v).The starting concentration for the oxidants ranged from 100 to 1000 µM, and every compound was measured in a series of measurements with 1.5× dilution for 6-10 different concentrations.Every measurement was performed thrice.

Presteady-State Kinetics
Pre-steady-state kinetics assays were performed under aerobic conditions using SX20 stopped-flow system (Applied Photophysics, Leatherhead, UK).The reduction of BsFNR by NADPH and its reoxidation by a quinone or oxygen was evaluated at 460 nm and 600 nm.The reaction intermediate spectra were recorded at different wavelengths in the 450-750 nm range over various timepoints.The BsFNR in syringe 1 (10 µM) was mixed with the contents of syringe 2 (100 µM of NADPH or 100 µM of NADPH and 66-500 µM of duroquinone).All measurements were performed in triplicate.

Photoreduction of BsFNR
The photoreduction of BsFNR (20 µM) was performed under anaerobic conditions in 0.02 M Hepes/NaOH buffer solution (pH 7.0) using 5-deazaFMN (2 µM) and EDTA (8 mM) as photosensitizer and sacrificial reagent, respectively.The reaction mixture contained 10 mM of glucose, 50 nM of glucose oxidase, and 50 nM of catalase.Before the introduction of the concentrated enzyme stock solution, the solution in the sealed spectrophotometer cuvette was flushed with oxygen-free argon for 30 min.The 5-deazaFMN solution and cuvette contents were protected from light before illumination.The cell was subsequently illuminated for one-minute-long intervals using a 100 W incandescent lamp at a distance of 20 cm.The progress of the reaction was followed spectrophotometrically for 1-3 h; each spectrum is an average of two.

Figure 1 .
Figure1.The 3D structure of a BsFNR dimer (PDB id: 3LZX).FAD-binding domain is colored yellow, while the NADP(H)-binding domain is colored red.The box denotes amino acids comprising each domain in a monomer.The bound FAD is shown in purple.The molecular graphics were made using UCSF ChimeraX[14].

Figure 2 .
Figure 2. The determination of the redox potentials of BsFNR.The photoreduction of 20 µM of BsFNR at different times of illumination, with the arrows indicating the changes in the spectra.The numbers next to the curves correspond to the total duration of illumination in minutes.Inset shows the interdependence of absorbance changes at 457 and 600 nm during the photoreduction of 20 µM of BsFNR.

Figure 2 .
Figure 2. The determination of the redox potentials of BsFNR.The photoreduction of 20 µM of BsFNR at different times of illumination, with the arrows indicating the changes in the spectra.The numbers next to the curves correspond to the total duration of illumination in minutes.Inset shows the interdependence of absorbance changes at 457 and 600 nm during the photoreduction of 20 µM of BsFNR.

Figure 2 .
Figure 2. The determination of the redox potentials of BsFNR.The photoreductio BsFNR at different times of illumination, with the arrows indicating the changes in t numbers next to the curves correspond to the total duration of illumination in minut the interdependence of absorbance changes at 457 and 600 nm during the photoredu of BsFNR.
phenanthrenequinone, the calculated single-electron flux was equal to 70%, with superoxide dismutase (SOD) (100 U/mL) inhibiting cytochrome c reduction by 43%.The assessment of the single-electron flux for the reduction of ArNO 2 is also based on the ArNO 2•− -mediated reduction of added cytochrome c.In the presence of 50 µM of NADPH and 100 µM of TNT or p-nitrobenzaldehyde, the rate of BsFNR-catalyzed reduction of 50 µM of cytochrome c was equal to 180-200% of the NADPH oxidation rate, corresponding to 90-100% singleelectron flux.SOD inhibited the reactions by 10-30%, reflecting the rapid reoxidation of the radicals with O 2 and the participation of the superoxide in the reduction of cytochrome c.Int.J. Mol.Sci.2024, 25, x 7 of 17

Figure 5 .
Figure 5. Dependence of the reactivity of quinones, nitroaromatic compounds, aromatic N-oxides, and benzyl viologen on their single-electron reduction midpoint potentials (log10 scale).Numbers and reduction potentials of the compounds are given in Table1.

Figure 8 .
Figure 8.The kinetics of BsFNR (5 µM) reduction and reoxidation under multiple turnover conditions with a varying concentration of duroquinone in the presence of 50 µM of NADPH followed at 460 nm (solid curves) and 600 nm (dashed curves): 1-250 µM of duroquinone, 2-166 µM of duroquinone, 3-111 µM of duroquinone, 4-74 µM of duroquinone, 5-49 µM of duroquinone, and 6-33 µM of duroquinone (all concentrations were reported after mixing).The inset shows the Lineweaver-Burk plot for the dependence of the apparent first-order rate constant of the reoxidation of reduced FAD on the concentration of duroquinone.

Figure 8 .
Figure 8.The kinetics of BsFNR (5 µM) reduction and reoxidation under multiple turnover conditions with a varying concentration of duroquinone in the presence of 50 µM of NADPH followed at 460 nm

Figure 9 .
Figure 9.The spectra of reaction intermediates formed during the turnover of 5 µM of BsFNR in the presence of 50 µM of NADPH and 250 µM of duroquinone (concentrations were reported after mixing).Differences in absorbance are shown at several timepoints over the 450-750 nm range.